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Abstract 

We present an entirely microscopic calculation of the Casimir force 
f{d) between two metallic plates in the limit of large separation d. 
The models of metals consist of mobile quantum charges in thermal 
equilibrium with the photon field at positive temperature T. Fluc- 
tuations of all degrees of freedom, matter and field, are treated ac- 
cording to the principles of quantum electrodynamics and statistical 
physics without recourse to approximations or intermediate assump- 
tions. Our main result is the correctness of the asymptotic universal 
formula f{d) ~ — ''^g^^f'^ , d — > cxd. This supports the fact that, in the 
framework of Lifshitz' theory of electromagnetic fluctuations, trans- 
verse electric modes do not contribute in this regime. Moreover the 
microscopic origin of universality is seen to rely on perfect screening 
sum rules that hold in great generality for conducting media. 

1 Introduction 

Motivations 

In 1948, Casimir [1] predicted that two neutral metallic plates placed 
in vacuum at distance d attract one another due to the electromag- 
netic field's zero-point fluctuations. In his calculation the microscopic 
structure of the conductors is not taken into account. The latter are 
treated as macroscopic bodies imposing metallic boundary conditions 
on the Maxwell fields. 
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Early experiments aimed at probing this theoretical prediction re- 
mained inconclusive [2] until the late 1990s where first experimental 
demonstrations were performed [SI Ul |5], opening the way to many 
others (see, e.g., |6l Sec. 3.5] and [7] for short reviews). A quanti- 
tative comparison with experiments requires to include a number of 
effects not accounted for in Casimir's simple treatment, such as the 
finite conductivity of the plates, the roughness of the surfaces and the 
dependence on the temperature T. 

Lifshitz [SI |9] provided a first major generalisation by consider- 
ing plates whose electric properties are described by a frequency- 
dependent dielectric function e{uj). The plates are in thermal equi- 
librium with a stochastic electromagnetic field, whose random nature 
is generated by the quantum and thermal fluctuations of photons and 
matter. The general force formula obtained in this way covers in prin- 
ciple a broad diversity of media, and should be valid for all regimes of 
temperature T and plate separation d. These regimes are character- 
ized by the single dimensionless parameter 

he Aph , , 

which measures the ratio of the photon thermal wavelength Aph = Phc 
to the separation distance d {h is the Planck constant, c the speed of 
light, fee the Boltzmann constant, /? = l/k^T). 

Although the Lifshitz theory (along with its various reformulations 
|1U| [TT1 [T21 El E!) is commonly used to interpret the experimental data, 
its predictions are uncertain when applied to conducting media at 
nonzero temperature. Indeed, the value of the force then depends cru- 
cially on the behaviour of the dielectric function at vanishing frequen- 
cies — a behaviour not directly accessible to experiments. This has 
led to several theoretical and numerical studies, resulting in a debate 
that has not yet evolved to a concensual end (see [3 |T3l HH |T5l [16] and 
references cited therein). In short, the controversy amounts to know- 
ing whether the reflection coefficient of the transverse electric (TE) 
mode of the field, r™(u;,fc) (depending on e(w)), vanishes or not in 
the limit of zero frequency. In the low-temperature, small-separation 
regime a ^ 1 one finds the force by unit surfac^ 

/(<i)~-|i^ + 0(T*), ifrTE(0.fc) = l, (2) 



^In these formulas, negative/positive terms stand for attractive/repulsive contributions. 
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while in the high-temperature, large-separation regime q <C 1, one 
has: 



/M--^!^, ifrTE(0,fc) = 0. (5) 

In Q-Q, the dominant term is in both cases the standard Casimir 
result, whereas in the high-temperature, large-separation regime Q- 
([5]), one sees a striking reduction of the force amplitude by a factor 
1/2 when the TE field modes are assumed not to contribute. Let us 
add that Formulae Q can be obtained under the assumption of 
the plasma relation for the low-frequency dielectric function (e(a;) ~ 
1 — cjp/u;^, u; — > 0, cOp the plasmon frequency), whereas the Drude 
expression (e(a;) ~ l—Airia/uj, a the conductivity) leads to the results 
([3]), ([5]). The force Q is also retrieved by extending Casimir 's original 
calculation to finite temperature between macroscopic plates that are 
not subject to charge fluctuation [T71 [T8] . 

Most of the actual debate focused on the finite-temperature correc- 
tions in the low-temperature, small-separation regime. In particular, 
the fact that the corrective term linear in T in ([s]) reflects a nonzero 
entropy at T = consists in an unacceptable violation of the Nernst 
postulate for some authors [Tj [H] , who thereby favor a nonvanishing 
reflection coefficient r™(0,fc). Other authors |15tll9]. however, argue 
that this linear correction no longer holds at very low temperature, 
and favor 

In this paper, we pronounce on the controversy in the large-separation 
regime (with fixed positive temperature). In order to decide which of 
the two alternatives Q and ^ is correct, we present a fully micro- 
scopic treatment of the Casimir effect based on the principles of quan- 
tum electrodynamics and statistical mechanics which does not suffer 
of intermediate models, assumptions or approximations]^ By fully mi- 
croscopic treatment we mean that all degrees of freedom, matter and 
field, are taken into account, contrary to Casimir 's original calculation 
that ignores particle fluctuations inside the plates. Microscopic mod- 
els have been produced to retrieve and justify Lifshitz' formula in the 
case of dielectric matter [231 IMl ES] j but conducting media offer more 
difficulties as far as one has to deal with screening phenomena due to 
free charges and magnetic forces between free currents. 



^Experimental setups in cylinder-plane and parallel plate geometries are currently being 
developed with the purposes of discriminating between the different proposed values for 



the force [Mllllllll]. 
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In [26], we computed the average force by unit surface between 
slabs containing purely classical charges and interacting via the static 
Coulomb potential, finding 

/('^) "^^^^ (6) 

In the letter [27] we improved the calculation by considering slabs 
made of quantum charges and interacting, in addition to the Coulomb 
force, with the transverse part of a classical electromagnetic field. 
These features do not alter the form ^ of the large-separation asymp- 
totic force. Reference [27] together with the companion letter |28j 
stress the importance of including in the calculation the effects of the 
charge fluctuations in the metals, which are responsible for reducing 
the asymptotic force amplitude Q by the factor 1/2. It is also un- 
derstood why the entirely classical model treated in [26] correctly pre- 
dicts the high-temperature, large-separation result (|6]): this is a con- 
sequence of the Bohr-van Leeuwen theorem. The theorem states that 
in classical systems at thermal equilibrium, matter decouples from the 
transverse electromagnetic field. Since high-temperature conducting 
phases tend to behave classically, the corresponding Casimir force will 
be determined at leading order by purely electrostatic interactions. 

Statement of results 

The present paper addresses the question of the Casimir force in the 
general framework of nonrelativistic thermal quantum electrodynam- 
ics (TQED), namely nonrelativistic quantum charges in interaction 
with the quantized electromagnetic field. The model (described in 
more detail in Section 2) consists of mobile quantum charges con- 
fined in two slabs A and B of thickness a and b with lateral faces of 
surface L^, set at distance d from each other. The charges interact 
with a quantum electromagnetic field enclosed in a large box A. The 



Hamiltonian -ffA,L,d of the system is specified in Formula (10) of Sec- 
tion 2. The photons and the particles are supposed to be in thermal 
equilibrium at temperature T, so that all the relevant information is 
contained in the grand-canonical potential $A,L,d associated with the 
Hamiltonian Hj^ j^ d. The average force by unit surface exerted be- 
tween the plates is defined by the rate of change occasioned in ^>A,L,d 
when varying the separating distance d: 

1 d 

The Casimir force is defined as 
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where the thermodynamic hmit of the system is taken in two stages. 
We first let the box enclosing the field A ^ M^, and then extend the 
plates' surfaces L? — > M^. The plates' thicknesses a and h are kept 
finite. The main result is that at any fixed temperature T > (such 
that the thermal energy ksT is much less than the rest mass energies 
rrK? of the particles), the large-separation asymptotic force is again 
given by Formula (|6]). The amplitude is linear in T, independent of 
the Planck constant h and of the speed of light c, and universal with 
respect to the microscopic constitution of the plates. Nonuniversal 
contributions and contributions depending on h and c, will only oc- 
cur at the next order 0((i~^) in the large-separation expansion of the 
force. This result not only validates from first principles the second 
alternative in Eqs. Q, ([5]) (associated with the vanishing of the re- 
flection coefficient for TE modes), but also establishes universality of 
the leading term ([g]) on a microscopic basis. Since our methods might 
be not quite familiar, we summarize the main steps of our derivation. 



Methods 

Casting the quantum system in a classical-like form: the 
space of loops 

When we have a conducting medium it is of utmost importance to 
deal properly with the collective screening effects. The idea is to 
cast the quantum system in a form which is as close as possible to 
that of a classical system of charges, where we have for instance the 
well-developed Debye-Hiickel theory of screening. To this effect we in- 
troduce a joint functional integral representation of the Gibbs weight 
associated with the total Hamiltonian of matter and field. In this for- 
malism, developed in [29] and recalled in Section 2, quantum particles 
appear as extended objects consisting of random closed wires C (called 
loops) carrying both a charge and a current. The loop size, which is 
measured by the thermal de Broglie wavelength, refiects the intrinsic 
quantum fiuctuation of the particle. There are two kinds of pairwise 
interactions between loops. The first one V^{Ci, Cj) originates from 



the electrostatic (Coulomb) potential between charges. Formula (20|. 
The second one, W™(£i,£j) (Formula (|2T|)) is called the magnetic 
potential. It is an effective interaction resulting from integrating out 
the field degrees of freedom: one can figure it as current interactions 
between the loops mediated by the transverse part of the electromag- 
netic field. At this point, although being an exact representation of 
the quantum TQED system, the statistical mechanics of loops has a 
classical-like structure which enables a convenient application of the 
methods of classical statistical mechanics. 
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Expressing the Casimir force in terms of loops 

The force between two loops is given, as in classical physics, by the 
gradients of the potentials dxV^{Ci, Cj) and dxy\^™{Ci, Cj) (along the 
X-axis perpendicular to the plates). The average force, as usual, is 
obtained by averaging these forces with the equilibrium correlation 
function p^'^\Ci, Cj) between two loops. The precise expression is 



found in Formula (38) in Section 3 where we have also singled out 
the proper Casimir force due to fluctuations. The additional part, 
called here capacitor force, is the direct Coulomb force that would 
occur between globally nonneutral plates. One can benefit from the 
translational invariance along the y directions parallel to the plates by 
using the correponding two-dimensional Fourier variable k and scale 
it as k = q/d where q is now a dimensionless Fourier variable. The 



scaling trivially provides a prefactor l/d? in the force Formula (40). 
The remaining d-dependence remains embedded in the microscopic 
expressions of the forces and correlations between loops. 

Screening of the electrostatic interactions 

This is the subject of Section 4. The main observation is that the 
Coulomb interaction between loops can be decomposed into = 



yei _^ (ggg Formulas (42 ) and (43 )). Here is the genuine clas- 



sical electrostatic interaction between charged wires whereas W"^ in- 
corporates the proper effect of the quantum nature of the particles 
manifested by the fluctuations of the loops. One can easily extend 
to the standard ideas of the classical Debye-Hiickel theory, pro- 
viding an effective resummed potential $ that becomes integrable at 
large distances (in the planar geometry one is concerned by the inte- 
gr ability in the y directions along the plates, see Appendix C). One 
is therefore left with the screened potential $ together with the addi- 
tional interactions VV"^ and W™. The latter would not be present for 
classical charges: they result from the intrinsic quantum fluctuations 
of the particles and behave as electric and magnetic dipole interaction 
at large distance. At this point on can use in the space of loops the 
methods of Mayer expansion and integral equations well developed in 
the context of classical Coulomb fluids. Of particular importance is 
the perfect screening sum rule stating that any specified loop is sur- 
rounded by a screening cloud of loops whose total charge compensates 
that of the specified loop. This imposes an exact integral constraint 



on the two-loop correlation function. Formula (52), that turns out to 



be at the origin of the universality of the Casimir force. 
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The large-separation asymptotics 

Apart from the obvious l/d^ scaling factor, the d-dependence of the 
force has to be extracted from the forces and correlations between 
loops. The Coulomb part of the force (expressed in terms of the 
q transverse Fourier variable) has a non vanishing limit as d ^ oo 
whereas the magnetic force vanishes as and the correlation be- 
tween the plates as 1/d (Section 6). This implies that the dominant 
term decays as 1/d? and does not involve a direct contribution of the 
magnetic force. To obtain the amplitude of this 1 /d? term it is neces- 
sary to determine the exact asymptotic form of the correlation. The 



latter is seen to be made of two terms (Formulae ( 59 ) and ( |62[ ) ) de- 
caying as 1 /d times a product of certain arrangements of correlations 
pertaining to the individual plates. When this is introduced in the 



force Formula ( 40 ) one discovers that the perfect screening sum rules 
in each of the plates wash out all details of the microscopic structure 
of the conductors, thereby leading to the wonderfully-simple result 
([g]) and providing a physical explanation of universality (Section 5). 
The analysis of the correlation in Subsections 6.1 and 6.3, if somewhat 
lengthy, uses common reasoning in terms of Mayer graph diagramatics. 
It basically reveals that the electric and magnetic dipolar potentials 
between the loops W'^ and W™ do not eventually contribute to the 
dominant term of the force. This term is entirely due to elec- 
trostatics and screening, so explaining why the purely classical model 
of [26] gives the correct result. The complexity of the full quantum 
mechanical treatment presented in this paper contrasts with the sim- 
plicity of the result, a fact that we could not foresee right away. More 
comments and perspectives are offered in the Concluding remarks in 
Section 7. 



2 Description of the system 

We consider two parallel slabs A = [—a, 0] x and Bd = [d, h+d\ x 
with thickness a and h and lateral surface L?. The x-axis is perpen- 
dicular to the plates, the inner face of slab A being fixed at x = 
while the slab is set at a distance d from it. The slabs contain 
nonrelativistic point particles of several species 7 (electrons, ions, nu- 
clei) with charges e^, masses m-y, spins and appropriate statistics. 
These particles are confined by walls without electrical properties in 
the two separate regions and no exchange is possible from one slab to 
the other. Particles in one plate are always distinguishable from alike 
particles in the other plate. To ensure the global neutrality of each 
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plate, we impose 



a b 

where the sums are carried over particles in A and i?^, respectively. 

This system of interacting charges is coupled to a quantum elec- 
tromagnetic field which is itself enclosed into a larger box A englobing 
both plates. The A^-particle Hamiltonian reads 

^1 e 2 

4=1 i<j 
N 

i=l 

with v{ri — Tj) the static Coulomb potential 

v{ri-rj)= ^ (11) 

As is common in atomic physics when matter is nonrelativistic and 
high-energy processes are neglected, we use the Coulomb gauge and 
electrostatic Gaussian units |30j. The Coulomb gauge has the advan- 
tage of clearly disentangling electrostatic and magnetic couplings in 
the Hamiltonian. The divergence-free vector potential A{r) is sup- 
posed to satisfy periodic boundary conditions on the sides of the box 
A. Its expansion in Fourier modes K is given by 

A(r) = (^)'^'j:5(^)^(ak,Ae-^^-'- + aK,Ae^^-n> (12) 

where ^ , ax, A are the creation and annihilation operators for the 
mode K, A of frequency ujk = c\k\ with commutation relations [aK,X7 x' 

k:"^a,a'; ^k,\-, a = 1, 2 are the polarization vectors; ^(a"), ^(0) = 1, 
is a real, spherically-symmetric, and smooth form factor taking care of 
ultraviolet divergencies. It is supposed to decay rapidly to beyond 
the characteristic wavenumber Xcut = — jt. where fh is an average 
particle mass. The term m""^ in (|10|) is the free field Hamiltonian 



^0A=E^^^«*^,A«K.^- (13) 

The wall potential 7i) confines the particles either to slab 

A or to slab Bd-, depending on whether 7^ designates a species in A 
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or Bd- Note that we neglect spin-field couplings in this model (see 
comments in the Concluding remarks). 

The states of this system of particles and field are supposed to be 
thermalised at the inverse temperature (3 = {kBT)~^, and statistical 
averages, denoted by (. . .), are taken with the usual Gibbs weight 
q-I3Ha,l_ ■\Ye introduce the finite- volume grand-canonical potential of 
the full system 

$A,L,d = -feBTlnTYe-^(^A.i-'^-^) (14) 

where the trace Tr = TrmatTrrad is carried over particles' and field's 
degrees of freedom. Here /x = {/^7„,^7i,} is the collection of chem- 
ical potentials that fix the average particle densities in each of the 
plates and AT" = {N^a^-^-yb} corresponding particle numbers. 



In (14) Trmat is carried out only on neutral configurations in each 
plate. The average force by unit surface exerted between infinitely 
extended plates immersed in the electromagnetic field is then defined 
by Formula ([s]), the temperature and chemical potentials being fixed. 

In addition to the separation d, there is a number of other char- 
acteristic lengths in the system, in particular, the thermal wavelength 
of photon Aph = /She and of particles Amat = hy^(3/m. Moreover, 
the plates are assumed to be conducting, and therefore characterized 
by a screening length Ascreen- Our derivation holds for the following 
hierarchy of lengths: 

Acut = ^^"^ ^ < Amat < Aph = V /3mc2 Amat < d, (15) 



mc 



Ascreen <. a,b <^ d. (16) 



The first set of inequalities is necessary for the consistency of the 
nonrelativistic treatment of matter, which requires (3fhc^ > 1 (the 
thermal energy is much smaller than the rest mass energy of the par- 



ticles). Inequality (16) means that the plates' thickness should be 
large enough for allowing the screening mechanisms to take place in- 
side the conductors. Finally, the conductors will be assumed to be 
invariant under translations and rotations in the plate directions. 



Loop formalism 

Our analysis relies on the formalism developed in [29] , based on a joint 
functional representation of both matter and field. In this formalism, 
the field degrees of freedom can be integrated out exactly. Then, the 
particle variables live in an auxiliary classical-like phase space whose 
elements are loops of random shape (for the statistical mechanics of 
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charged loops, see the review [31], Chap. V and references therein). 
A loop £ = (r, x) is specified by a position r in space and a number 
of internal degrees of freedom x = {l,?, ^i')) consisting of a species 
index 7, a charge number p £ 1,2,3, ... and a closed Brownian path 
s I— > X{s), s £ [0,p], -X^(O) = X{p). The loop's shape X(s) is a Gaus- 
sian stochastic process (Brownian bridge) whose functional integral 
has unit normalization, zero mean, and covariance given by 

jiy{X)X^{s)X%s') = V p (min ^} - . (17) 

The loop's path is 

rW=r + A^X(s), 0<s<p (18) 

where = h^J~^JrrUj is the de Broglie thermal wavelength. The oc- 
currence of the Brownian path results from the Feynman-Kac path 
integral and gives the extension of the quantum particle's fluctua- 
tion. The number p accounts for the quantum statistics of the species 
7. It corresponds to grouping together p particles that are permuted 
accordingly to a cyclic permutation of length p. 

The pairwis6 interaction 6'y^6'y^ V^(i, jf) between two loops % — Ci 
and J = Lj is the sum of two contributions 

e,,e^^.y(i,j) = e^,e^Jy=(i,j)+W'"(i,j)]. (19) 

The first contribution, inherited from the Coulomb potential, is 



V'{i,j) = jisijis,5{si-s,) (20) 



where s = s mod 1 and 5{si — sj) takes into account the equal-time 
constraint imposed by the Feynman-Kac formula. The second contri- 
bution is the effective potential resulting from the elimination of the 
field's degrees of freedom. We call it the magnetic potential. It is 
given in Fourier representation by Formula (66)-(67) of [29]: 

^"^^'■^'^ = IW?^'""^^^'^'^ W-(xi,X,,i^), (21) 
W-(X^,X,,K) = ^ rkxf(..)e^^-^-^»(-) 



X 



rdXj(.,) e-^-H (^.) '-^CiK) QiK, - s,), 
Jo ^ 



with 



Q{K,Si-sj) ^ . Z';!''^^ cosli[Xp^K{\si-Sj\ - 1/2)], (22) 
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and 



5'^,{K) ^ V - ^ = {K}l=i^ ^ = 1^1) (23) 



is the transverse Kronecker function. In (21 ), dX^{s) are stochastic 
line integrals along the loop shape. The function Q, depending only 
on Aph, is the manifestation of the quantum photon field. This formula 
holds when the field region A has been extended to infinity, replacing 
the discrete sum on Fourier modes by an integral. 

Written in terms of loop variables, the grand-canonical partition 
function of the full system, normalised by that of the free radiation 
field, has a classical structure: 

rp^ g-/3(/fA,i,-M-Ar) 
S r ,y = lim 3 

oo oo » n 4 „ ng 

nA=0 ^■n£j=0 ^'-^^ a ''^d 

(24) 



In ([24]), the loop integration f^dCa = /^dr^ dxa = Xl'i'^a Y.-yaY.pPi^a) 

is carried over paths rjj entirely contained in slab A, respectively, 
J^^dCb over paths in slab B^i. This corresponds to choosing hard 
walls on the faces of the slabs, i.e., Dirichlet boundary conditions 



for the particle wavefunctions. In (24) the sums run only on neutral 
configurations of loops in each slab. 

The total loop energy U can be separated into intra and interplate 
contributions: 



U = Ua + Ub, + Uab, (25) 



where 



Ua=J2 e7,e7.n^>i) (26) 

is the sum of interactions occurring among loops confined into slab A 
(likewise for Ub^ in slab B^), and 

Uab, = Y.H ^y.^^.mJ) (27) 

is the interaction energy between the two plates. Moreover, each loop 
is equipped with an effective activity z{jO) containing the loop self- 
energy e2y(£,£): 

p (2vrpA2)3/2^ • ^^^> 
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The factor (2s^ + 1) accounts for the spin-degeneracy of the energy 
levels, /i-y is the chemical potential of species 7 and r/^ = ±1 for 
bosonic/fermionic species. 

We stress that although the loop partition function ( 24 1 has a clas- 
sical form, it is a mathematically exact representation of the original 
grand-canonical partition function of the system of quantum charges 
and photons as defined by the Hamiltonian (10). 



3 The Casimir force 



In view of (14)-(|8|) and the fact that the free-field partition function 



does not depend on d, the Casimir force expressed with the help of 



the loop partition function ( 24 ) reads 
kuT d 



f{d) = lim lim 



InTr e-/5(^A,^-/x-iV) _ i^TY e'^^S^A 



L^ooA^RS L"^ dd 
lim ^ (29) 



The dependence upon d in the partition function (24) occurs only in 
the confinement of the loops in the slab B^. For the rest of the paper, 
it is convenient to shift the positional integration variable S [d, b+d] 
of a loop in to Xf, — d, so that slab B^ is moved to the fixed region 
B = [0, b] X L^. Then, the d-dependence is transfered to the interaction 
potential between slabs: 

V{Ca,Cb), xbe[d,b + d] ^ V{Ca,Cb + d), xbe[0,b], (30) 

where £f, -|- d is the loop £ft shifted along the x-axis from to Xh + d. 
This amounts to measure the positions in slab B^ from its inner face. 
To abbreviate the notation, we set 

VAB{Ca,Cb) = V{Ca,Ct + d). (31) 

From now on it will be implicitly understood that VAsi^a, ^b) depends 



on d according to (31) and that in forthcoming integrals the path of 
the loop Ca is restricted to the fixed slab A of volume [—a, 0] x 
while that of Ct to the fixed slab B of volume [0,6] x L^. 



Differentiating with respect to d in ( 29 ) is equivalent to differenti- 
ating the potential V{Ca,Cb + d) with respect to Xb, or with respect 
to —Xa (since the dependence on the x components is Xa — Xh — d). 
This brings in the average force along x 

f[d) = ^im -1 / ^e^,e^,(a,„yiB + d,^W^B){C-a. cJS • 

\ a b I loops 

(32) 
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The bracket (• • •)ioopg denotes the grand-canonical statistical average 



in the phase space of loops with activities (28) and with respect to 
the Gibbs weight e~^^ associated to the loop potential energy (25). 
Since the force is a two-body observable, its thermal average can be 

(2) 

expressed as an integral over the two-loop correlation p)^ between a 
loop in A and a loop in 

fid) = hm j^[dl[d2e,,e,,{d,,VXB + d.,W^B)i^,2)p'-2lL{l,2), 

(33) 



where we have followed the notation (31 ), 

p(,%^(£i,£2)^pf (/:i,/:2 + d), (34) 

and the dl integration is carried on loops in A, respectively, the d2 
integration on loops in B. 

At this stage we take the limit L — > oo of infinite plate surfaces. 
Since the conductors are assumed to become homogeneous in the y = 
(y, z) plane of the plates, the two-loop correlation function tends to 

(2) 

a function ^0)^^(1,2, lUi — y2\)- In this geometry, it is convenient to 
decompose ri = {xi,yi), and 

l = (l,yi), dl = dldyi, dl = d2;idxi (35) 

where 1 = {xi,xi) denotes the position along x of the loop £i and its 
internal degrees of freedom (likewise for 2). In the limit, the factor 
1 /L^ cancels with one of the ?/-integral in ( 33 ) , yielding 

m = /dl/d2/dy e^,e^,(5,.,y|B + 5,,W^B)(l,2,y)pj;j(l,2,y). 
J A JB J 

(36) 



We introduce in (36) the loop Ursell function h{l,2) defined in the 
usual way by 

p{l)p{2)h{l, 2) ^ 2) - K1)P(2), (37) 

with p(>C) the loop density, so that 

/(d) = /dl/d2/d?/ (5,,yi5 + 9,,WXB)(l,2,2/) 

J A JB J 

X /5A(l)pB(2)/iAB(l,2,y) + /cap(d), (38) 



3 



Loop correlation functions are defined similarly to particle density correlation func- 
tions, see [Slj, Chap. V. 
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where again hAB{^i,^2) = h{Ci,C2+d) andp^(£i) = p(£i), ^^(£2) 



p{C2 + d) according to the notation (31), (34). The capacitor force 
fca.p{d) = /cap + fclp comes from the subtracted product of the plates' 
density in (37). We show in Appendix [A| that the electrostatic part 
/cap(^) due to the term S^i^ab reduces to 



/cap(rf) = 2vr 



dxi ca{xi) 



dX2 Cb{x2) 



(39) 



where ca{xi) and cb{x2) are the mean charge densities in plate A and 
B. It corresponds to the standard force (in Gaussian units) between 
a capacitor's plates whose surface charge densities are respectively 
/°^dxiCyi(xi) and jQdx2CB{x2)- In this work, we assume strict neu- 



trality in virtue of (|9|), so that fcapi^^) = 0- Moreover, we also show in 



Appendix [A| that the magnetic contribution f^p{d) decays faster than 
any inverse power of d. In the sequel, we will thus drop the capacitor 



force and focus only on the first term of (38), which is the proper 
Casimir force generated by fluctuations. 

There is a noteworthy simplification in the electrostatic part due 
to 9x1^^(1) 2, y) in the Casimir force. Namely, one can omit all mul- 
tipolar contributions in the loop Coulomb force c^xiC-j-j^e-^j 

y=(l,2) (see 



(20)) replacing it by its pure monopole term 9x1^1^2671672 |ri-r2| ~ 
PiP2e'yie-y2dxiv{ri — r2). The underlying reason is that the electric 
force expressed in terms of the original two-point particle correlation 
function involves the average standard Coulomb force dxiv{ri — r2) 
between point charges. The equivalence with the present formulation 
in terms of loops is given in Appendix [B] 

Finally, we represent the remaining ^/-integral in the two-dimensional 
(transverse) Fourier space k and introduce the dimensionless variable 
q = kd: 



fid) 



1 

1^ 



dl/d2 

A Jb 



dq 



(2vr)2 

XpA{l)pB{2)hAB{l,2,^), 



671 672 {PlP2 dx^VAB + Qx^Wab) (1> 2, 2) 



where 

5x1^^5(1,2, f) 



e'3 



V^(xi - X2- dY + y2 



(40) 



27re"'?e"^(''2-xi)/d_ 



(41) 



The general formula (40) is an expression structurally similar to the 



one developed in the purely classical model |26| Form. (29)] . It reduces 
to it when charges are classical and the field is switched off. The main 
purpose is now to extract the d-dependence of the Ursell function 
HabO-, 2, 3), which embodies all correlations between the two plates. 
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4 Screening of the electrostatic inter- 
action 



The Ursell function can be conveniently analysed by performing a 
Mayer expansion in the phase space of loops. The Mayer bonds 
f{i,j) = — 1 are built from the basic loop-loop interac- 



tion ( 19 ) . The Coulombic part of ( 19 ) decays as r -"^ so that the Mayer 



bond is not integrable. In order to remedy to this nonintegrability, it 
is necessary to take screening effects into account. To this end, we first 
make the following observation. From the Feynman-Kac formula the 



potential (20) inherits the quantum-mechanical equal-time constraint: 
i.e., every element of charge e^- A^-dXj(si) of the first loop does not 
interact with every other element e-y^X^^dX j(sj) as would be the case 
in classical physics, but the interaction takes place only if si = S2. It 
is therefore of interest to split into V^^ + W^, where 

V%,J) = Jl,ll,-j;^ (42) 

r « j 1 

W'{i,j) = dsjdsj{6{s,-s,)-l) — ^ (43) 

Jo Jo \rr^ - ry^\ 



The contribution is the genuine classical Coulomb interaction be- 
tween two uniformly charged wires of shapes rf'^ and f^j^\ whereas 
the quantum-mechanical constraint appears in W^. 



Now, the complete two-loop potential (19) reads 

V = V^^ + W, with W = + W'". (44) 

It is known that and W™ have a dipolar r"^ decay at large dis- 
tance, see Section VI of [29]. One can therefore view the system of 
loops as behaving like classical random charged wires (interacting with 
V^^) with additional electric and magnetic multipolar interaction W. 

We deal with the screening effect generated by the classical Coulom- 
bic part V^^ by the standard Debye-Hiickel method. This amounts 
to introduce the effective screened potential $ corresponding to the 
chain-resummation of the linear part — /3e-y.e^^.y^(i, jf) of the bond 
f{i,j): ^ satisfies the integral equation 

= V%,j) - jdl'^V^, (45) 

where 

K-\l) = [4nl3elp{l)]-'/' (46) 
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defines a local screening length in the system of loops. This potential 
is now short-range in the sense that it is integrable on the y-direction 
along the plates (see Appendix C), implying 

lim |^>(l,2,fc)| < oo. (47) 

fc^O 

The Mayer series is reorganised by the Abbe-Meeron resummation 
process (|3I], Chap. V and references therein) into so-called "proto- 
type" graphs n with integrable bonds F{i,j) and F^(i,jf)|^ given by 



F{i,j) = -Pe^,e^^'^{i,j), (48) 
F^{i,j) = e-/3^..^..(*+>V)(ij) _ ^ ^ pe,^e^^Hi,j). (49) 

The resummed Mayer graph series of the Ursell function reads 

h{h2) = Y,^j<^^p{^)--j^rnp{m) \[ (50) 
n {i,j}en 

where J- S {F, F^}. The diagrams 11 have two root points and m — 2 
internal circles (m=2,3,...), and a symmetry number ^n. In (50), the 
weights of the integrated points are the density, so that the graphs 
contain no articulation points. Prototype graphs are subject to an 
important rule: convolution chains of bonds F are forbidden to avoid 
double counting of the original Mayer graphs. 



Perfect screening sum rules 

On the microscopic level, the conducting behaviour of a system at 
equilibrium is characterized by the fulfilment of the "perfect screening 
sum rule" [32]: a fixed charge in the system is neutralized by the 
mean charge density surrounding it. This property is expressed by 
the following constraint on the two-particle Ursell function: 

5:/dn..p(.,.70M...T...7.) = -.,. (51) 

71 

It turns out that the same perfect screening sum rule holds in the 
auxiliary system of loops 

jdlpie^,p{l)h{l,2) = -p2e^,. (52) 

The interpretation is the same: the fixed loop 2 with charge ^2673 is 
surrounded by a screening cloud of loops with opposite total charge. 

''At large distance, F^{i,j) ~ —/3e^^e^^yV{i, j) ~ r"^ is at the border of integrability. 
Hence, some care has to be exercised as it is the case in dipole gases. 
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The sum rule (52) holds in great generality for infinitely extended 



conductors, in particular for slab geometries. A justification of this 
sum rule is easily given when the loop Ursell correlation h is replaced 



by the single bond F. The equation (45) written in the transverse 
Fourier space reads 



<I>(i, j, k) = V%, j, k) - jdl^V-\i, 1, fc)$(l, j, fc), 



(53) 



where from ( 42 ) 



V^\i,j,k) 



Jo 



(54) 



Here X(s) and Y(s) are the components of X{s) along x and in the 
y plane and a;!*! = x + A-yX(s) is the component along x of r^'^l (18). 
We divide both members of (53) by V'^^{i,j, k) and let k ^ 0. In view 



of the fact that ^{i,j,k) remains finite (see (|47|)), that 
limi^^QV'^^{i,l,k)/V'^^{i,j,k) = pi/pj, and from the definition (46) 
one obtains 



y"dlpi e^, p(l)F(l,j,fc = 0) 



(55) 



which is the same as ( 52 ) with F replacing h 



In fact, it can be shown that the general case (52 ) is a consequence 



of ( 55 ) , by using the same dressing argument as that presented in 



(63)-(68) of [26, Sec. 5]. 



5 Asymptotic Casimir force 



To analyse the asymptotic d-dependence of the force (40), we need 
to extract that of the electrostatic part dx^v^B and of the magnetic 
part Sj^jW™^ together with that of the Ursell correlation Kab- It is 
immediate from (41) that ^^^^^^^(l, 2, %) has the limit 



dx,VAB{l,2,l) 



27re'9 = 0(1) 



as d ^ oo. We will establish in Section |6] the following facts: 

5,,W^5(l,2,2) = 0(O, 
/iAB(l,2,i) = 0((i-i). 



(56) 



(57) 
(58) 



As a consequence, the average of the magnetic part dx^W™Q{l,2, |) 
does not contribute to the Casimir force at leading order since it is 
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0(d^^), whereas the electrostatic part of the force is 0{d^'^). To cal- 
culate the coefficient of this oc d^^ dominant contribution, the exact 
structure of H ab at 0{d~^) is needed. The latter is analysed in detail 
in Subsection 6.3 In short, both bonds Fab (|48]) and (49) are of 
order d~^ and the diagrams contributing to /i^b(1)2, at this order 
comprise only one of these AB links. Those having a single Fab bond 
sum up to the factorized expression 



1 q G°(1,0,0)G°,(0,2,0) 
/3(i47rsinhg 



where 



(59) 



G°i(l,0,O) 



a(1,0,0) 



F^{l,i,0) + 



{h^AriiM (60) 



comprises internal correlations occurring in slab A. The superscript 
"0" qualifies statistical-mechanical quantities characterizing the sys- 
tem governed by the same Hamiltonian (10) but where V""^^^^ confines 
particles in a sing le slab (either A or B)71?^{1, 0, 0) = h°^{l, 0,k = 0) 
is the Ursell correlation between a loop "1" in slab A and a classical 
charge located at its right border, denoted by the loop argument 







(x = 0, ao, p = l, X{ 



0) 



(61) 



The structure of is determined by the excluded convolution rule 
applied to Fab- The partial Ursell function (/i^)""^ occurring in the 



right hand side of (60) is defined in Subsection 6.3 see (84). The same 



notations and definitions apply to the plate B. 

The diagrams having a sing le F^^ bond sum up to the expression 



di p^Ai^ Idj p%U) [4(1,^,0) + ^ 



X {-Pe,^e,^)WABii,j, f ) /i^(j,2,0) + M 



(62) 



where again and are the Ursell functions of the single plate 
systems A and B. 

The rest of the analysis relies on the application of the perfect sum 



rule (52) for loops. Indeed, introducing the contribution (62) into the 



force ( 40 ) , one builds the integral 



dl pie^.p'^Ai 



h'Ai 



l,i,0) + 



<5(M 

pA» 



(63) 



which vanishes by ( 52 ) 
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Introducing now the contribution (59) into the force (40), we see 
that the integrals on the two slabs factorize as 



CO 

dq 







sinhq 



From ( 60 ) , we have 



dip, e,, p\{lf^ = km e.. .VD^ 



(64) 



(65) 



The first equality follows from the sum rule (55) for the F bond in 
the single plate A. The second equality is again a consequence of the 
perfect screening for loops (52). Perfect screening in plate B implies 
similar identities for the second bracket in (64). 

Noticing that the g-integral provides the constant C,{2>) /2, the Casimir 
force at large separation is 



C(3) 



(66) 



which is the main result of this paper. 



6 Asymptotic correlations between the 
two slabs 

To extract the asymptotic large-separation behaviour of the Ursell 
correlation /i^b(1,2, |), we select the class of prototype graphs that 
give the dominant contribution by analysing them one by one, as done 
in [Ml App. C]. 

It is important to distinguish situations where arguments Ci , Cj 
both lie in the same plate or in the two different plates. As done 



before (see (31), (34)) we index any quantity with arguments Ci G 
A = [0, a] xW and £j G S = [0, 6] x with an index AB. We 
introduce a similar notation for interactions and correlations internal 
to a given plate, using the index AA, BB when loops lie in the same 
slab, e.(?., 

FAA{Ci,Cj) = F{£i,Cj), FBB{Ci,Cj) = F{Ci + d,Cj + d), 
hAA{Ci,Cj) = h{Ci,Cj), hBsid, Cj) = Kd + d,Cj + d). (67) 
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In the limit d ^ oo, the plates will have no mutual interaction any- 
more: yli?-correlations are expected to vanish whereas AA and BB 
quantities will tend to those pertaining to the system constituted by 
a single plate. Using the superscript "0" to qualify the statistical 
mechanical description of the single plates, one will have in particular 



FAAii,j)^F^{i,j), 



FBBii,j)^FUi,j] 



(68) 



In the next subsections, we analyse in more detail the behaviour of 



/iab(1,2, |) occurring in the force formula (40) at large d. 



6.1 Large-distance behaviour of the screened 
potential ^ab 

The main fact to be established in this subsection is the factorization of 
the screened potential ^AB{i,ji 3)) ^-t order d^^, into two independent 
parts pertaining to the individual slabs A and B. We extend the 
arguments developed for a system of classical charges presented in 
Section 3.2.3 of Ref. [33j. One observes first that this factorization is 
already present in the bare Coulomb potential V^^{i, j, k) (54 ). Indeed 



taking into account the shift (30) of the loops' positions as well as 



xp' < < xp- 



for all Si , Sj , one can write 



ke' 



-kd 



27r 





kd 



Pj 

dsj e 



-ik-X 



li ^3 (■^.7) 27r - 
1.. c 



y|:4(i,o,fc) viB{^,j,k). 



-kyjcp 



27r 



(69) 



We have identified the first bracket to the Coulomb potential inside A 
between a loop i and a loop variable denoted 0, corresponding to an 
(arbitrary) classical charge situated on the inner side of the slab (see 
(61)), and likewise for the second bracket. We show below that the 
factorization extends to the screened potential as d ^ 00 m. the form 



,1 



ABVl^J 



' d) 



ci>Vi,o,o) ^'^(o,j,o). 



(70) 



d An sinh q 

where <I>^(i,0,O) is the screened potential at fc = inside the single 
plate A between a loop i and a classical charge at its right boundary, 
and likewise for $^(0, j, 0). 
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In the chain summation of V^^ bonds that constitutes ^ab, ob- 
tained by iterating (53), we keep only the dominant chains, which 
turn out to be of order d~^. We follow the steps performed in the 
classical situation in |33J. We split every integral on internal convo- 
lution points into an A and a B contribution. Using again the no- 



tation ([67|, i.e. specifying as Ffyi, j, fc), y|'g(z, j, fc), j, fc) 

according to the location of its argument in A or B, V^^ chains are ex- 
panded into chains made of V^j^, Fj'g and V^g bonds. One notes that 
VX\{i,j,k) = ViB{i,j,k), and that Vi^^i^j^k) = VfB{j,i,-k) by 
space inversion in the y plane. We call V^bC^, J, k) a traversing bond, 
and chains that link A with B, traversing chains. Clearly, traversing 
chains that contribute to ^ABihji f ) have necessarily an odd number 



of traversing bonds V^'^. Let be the sum of chains containing 

exactly 2n + 1 traversing bonds. The contribution is a sum of 



■ ■ ^AA * ^AA * ^AB * ^BB * ^BB ' 



convolution chains of the type 

Using the factorization (69) of V^^, one can resum on eiter side of 
convolution chains of and Fj'g into quantities $^1^(^,0, 

and $bb(0, j, ^aa differs from the screened potential ^aa by the 
omission of traversing chains starting in plate A and returning to it, 
which describe part of the electrical influence ol B on A (likewise for 
^bb)- But for large d, these traversing chains do not contribute any- 
more and by (68), <&aa(^,0, and $b_b(0, j, also tend as d 



oo 



to the screened loop potentials 0, 0) and $^(0, j, 0) of the single 
plates systems. Hence, 



+00 qe 



2TTd 



(71) 



is of order d ^ with a factorized coefficient 
The contribution ^^'^'^ 

the type 



is then formed by convolution chains of 

4a * Kb *Vbb---- Using m, ^ and 
resumming again the convolution chains of V^^^ and on either 

extremity into ^aa and ^bb, one obtains 



'AB 



'^fBihJ,1) 



qe 



3 



2TTd 



'^AA{i,0 
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d2 



dl 



Air 



4(2) 
47r 



yj^^(0,l,|)c^BB(l,0,2) 



<^aa{0,2,1)vU2,0,1) 



^BB(0,J,f) 



(72) 



By definition of ^aa and ^bb, these quantities satisfy the integral 
relation (53) relative to A and B with inverse screening lengths k\{1) 
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and in place of The brackets in (72) thus reduce to 



$AA(0,0,2)-y|i^(0,0,i) 



2TTd 



27rd 



+ 0(1), 
+ 0(1). 



(73) 
(74) 



On the right hand side of ( 73 ) and ( 74 ) , the dominant terms come from 
the Coulomb potentials (see (54)), while the estimates 0(1) reflect the 



fact that ^y{j^(i,j,k) and ^BB{hj,k) are bounded in k (see ( |47[ ) and 
Appendix C) . This yields 



2iTd 



(75) 



By induction on n, one easily sees that $2;;+'^(i , j, |) receives instead 
a prefactor ^-^e"^"^. Summing over n = 1,2,3,... gives the result 



(70). It is interesting to see that the screened electrostatic interaction 
between the plates, at order d~^, involves only particles located close 
to the inner faces of the slabs. 

For later convenience we write the equivalent factorized form of 
the bond F disl). 



FABiiJ, f) 



d^oo 1 q F0(i,0,O) F0(0,j,O) 



f3d 47r sinh q 



(76) 



where and C/Jq are two charges located at the inner boundary of 
the slabs. 



6.2 Large-distance behaviour of the dipolar po- 
tential Wab 

The partial Fourier transform yVA_B(l,2,fc) is related to the three- 
dimensional Fourier transform W(xi,X2,-K') (-K" = {ki,k)) of (44) by 

WAB{l,2,k) = l^e^'''(^^-^'-'''>W{xi,X2,K), (77) 

remembering that WABi-^ii ^2) = VV(£i,£2 + d). Changing ki 
qi/d and setting k = q/d, 



WAB(l,2,i) = ^ /^e^''^^e-'^^W(xi,X2,t,f), (78) 



which shows that Wab(1,2,|) = 0{d ^) provided the integral has 
a limit as d — > 00. The analysis of W{xi,X2, k) at small K has 
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been carried out in [29 . It was observed that the dipolar electric 
part W"^ was screened by thermalized photons at large distance. As 
a consequence }V{xi,X2, k) behaves as X2; -f*^) when K 

0, where VK™ is the magnetic potential corresponding to a classical 
electromagnetic field (i.e., setting Aph = in W™). The behaviour of 
W^{xi:X2, K^) itself was worked out in |34| Form. (25)]: 

w"^{xi,x2,K) ~ /^^^^^ ^ rdx^{s^)rdx^{s2) 

X[K-X^{SI)][K-X2{S2)]^6]:,{K), K^O. (79) 

It is analogous to dipolar magnetic interaction between two classical 
current loops of shape -X^i(-) and X2{-)- Since this interaction only de- 
pends on the unit vector k = k/k, W(x1) X2, ^) asymptotically 
independent of d, implying 

WAB(l,2,2) = 0(d-i). (80) 



An explicit expression for the asymptotic form of (78) can be found 
in Appendix |Dj 

The potentials and W™^ could as well be separately anal- 



ysed in the same way. One sees on (78) that 9a:^>V™g(l, 2, ^) has an 
additional factor, so that 

WMl,2,1) = 0{d-'), d,,WMh2,1) = 0{d-'). (81) 
6.3 The Ursell function at order 0{d~^) 

From ([76]), (80 ) and the definition of (49), one sees that FABihji 5) 
0{d~^) and F^^{i,j, ^) = 0{d^^). It is clear that the decay rate of 
a prototype graph will depend on the number of its traversing bonds. 
A rough counting gives d~'^^^ d~'^AB where uab is the number of Fab 
bonds and the number of F^^ bonds. A closer inspection shows 
that this decay can be even faster, at least as 

^~2I ^-UAB ^-n\g ^ d^oo, (82) 
where /, < / < nAB+n\^—l, depends on the topology of the specific 



diagram and / = if there is a single traversing bond. Formula (82) 
can be established repeating the same steps as Appendix C of [2 
(with integrals over loop degrees of freedom) . 



From ( 82 ) , the slowest decaying Mayer diagrams have either uab = 
l,n5g = or UAB = O,?!^^ = 1, namely only one Fab or one 
bond. In forming the complete correlation function of the two-slab 
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system, these bonds have to be dressed at their extremities by appro- 
priate internal correlations of the individual slabs in conformity with 
the diagrammatic rules. Thus, the complete expression of the Ursell 
function /iyiB(l,2, |) at order 0{d~ 



IS 



hAB 



Daa * Fab * Dbb + D^a * ^ab * ^bb- 



(83) 



The formation of the dressing function D differs from that of 
because of the excluded convolution rule in prototype graphs: no Faa 
or Fb b bond can be attached alone to the extremities of Fab whereas 
there is no such restriction for F^^. 

The dressing function {D^g) of F^^ consists of all possible 
AA {BB) internal graphs. According to the discussion at the begin- 
ning of the section (see ([68|)), it tends to the Ursell function h\ (hP^) °f 



the individual plate. The corresponding contribution to /ia_b(1,2, |) 
is thus given at 0{d~^) by (62). The delta terms in ( |62[ ) account for 
the situation where no bonds are attached to the extremities of F^^. 

To deal with the excluded convolution rule when forming the dress- 
ing function DaAi we introduce the function h'^{i,j) defined by the 
sum of all prototype graphs that do not begin nor end with a F link 
alone. Its relation to the Ursell function is 



h = F + K''^*F + F*h"-''*F + F*hJ''' + h'' 



(84) 



Then Daa = h^^ + Faa 



^Ta + ^Ipa (likewise for Dbb)- With 



the factorization (76) of the link Fab, one sees that {Daa * Fab 



-Cbb)(1)2, %) has the factorized form (59) with 



Pa 



+ 



+ f: 



* {h'Ar 



Fl. 



(85) 



Since the bond Fab is already 0{d ^), all other quantities have been 



evaluated for single plate systems according to (68). The final form 
( 60 ) of G\ follows by noticing that the first three terms of ( 84 ) build 



as given in ( 85 ) . 



7 Concluding remarks 

In this paper the large-separation asymptotics ([6]) of the Casimir force 
between two conducting plates has been derived exactly from the prin- 
ciples of quantum electrodynamics and statistical mechanics for any 
fixed nonzero temperature, taking all microscopic degrees of freedom 
of matter and field into account. This does not give a direct proof 
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that the TE mode reflexion coeflicient does not contribute at zero fre- 
quency, but a strong evidence for it. The derivation apphes to any 
model of conductor consisting of mobile quantum charges. The latter 
can be negative and positive charge carriers (like ions and anions in 
electrolytes), or, e.g., form the one-component electron gas in the jel- 
lium model of a metal, the central common point to all these systems 
being the screening mechanisms and the perfect screening sum rules. 

Let us note that there can be no contradiction between the be- 
haviour ([g]) and the Nernst heat theorem. As mentioned in the intro- 
duction, it has been argued, and controversially debated, that the use 
of the Drude expression of the dielectric function (yielding (|6|) was 
not consistent with the Nernst postulate which requires that the en- 
tropy of the total system vanishes at zero temperature. In our setting 
the question arises in different terms. The asymptotic formula ([6| is 
definitely true, whereas the Nernst theorem is the separate affirma- 
tion that the QED Hamiltonian (10) has a unique (or not extensively 
degenerate) ground state, a nontrivial and uncorrelated mathematical 
problem. 

A number of questions deserve further studies. We have disre- 
garded paramagnetic forces due to the Pauli coupling —ficr-B of elec- 
tronic and nuclear magnetic moments //<t to the magnetic field B. 
Preliminary investigations using spin coherent states functional inte- 
grals indicate that such interactions result in an additional effective 
dipolar potential which, as the orbital diamagnetic part W™, will not 
contribute to the asymptotic force. 

We have kept the thickness a and b of the plates finite while the 
separation d — > cxd. Then, because of perfect screening, the asymptotic 
force is universal as well as independent of a and b. This corresponds 
to the present experimental situation where only thin coats of metal 
of order of 50 nm are deposited on a substrate [35]. Compared to 
separations ranging from 0.5 to 3 fim, the regime is clearly a,b <^ d. 
The opposite situation of thick plates a,b ^ d, namely taking here 
a = 6 = oo at the very beginning, requires a careful analysis since 
the magnetic potential W™(1, 2, |) looses integrability as c? ^ cxd (see 
the factor exp{iqi [xi—X2]/d) in (78)). Then x-integrals have to be 
performed before taking the limit d — > oo which appears to lead to a 
modification of the d~^ coefficient with small (nonuniversal) terms of 
order 0((/3mc^)~^)) . 

Expression ( 66 ) is the first term of an expansion in inverse powers 
of d whose terms will be of the form An/d"^, n > 4. The amplitudes 
An{p,T) are no more universal. They will depend on the thermody- 
namic and geometric parameters of the plates (temperature, densities 
p, thickness a, b) as well as their microscopic characteristics (particle 
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masses and charges). Looking at the form of the electrostatic and mag- 
netic dipole potentials W™ and W"^, the expansion can be cast in terms 
of dimensionless quantities including [Xm&t / d)"^ ■, (Aph/c^)", where Amat 
and Aph are the matter and photon thermal lengths. This expansion is 
therefore only meaningful when the condition (15) is met. Of great in- 
terest would be the calculation of the first subdominant amplitude ^4 
that includes corrections from an imperfectly-conducting metal and 
to compare it with the predictions of the Lifshitz theory [6j . Also the 
effect of the capacitor force, which cannot be completely turned off 
in experiment, can be estimated by analysing the term (39) at large 
separation. 

Finally, an open question is the understanding of the crossover 
to the zero-temperature Casimir force / ~ — vr '^hc/2AM'^ due to pure 
quantum fluctuations. In the framework of Lifshitz theory, it was 
shown that plasmon modes at the surfaces of the plates combine with 
photonic modes to build the above usual zero-temperature Casimir 
force calculated as if the plates were inert |36j . Notice that we have 
not added to the Hamiltonian (10) the vacuum energy ^'YliK \ ^K- 
In fact this (infinite) constant plays no role since it will anyway not 
appear in the force formula ([8| (it is independent of d in our setting) . 
To study the zero-temperature case one cannot rely on the above- 
mentioned expansion since Aj^iat ? Apji — > 00 as T — > and condition 
(15) does not hold anymore. One has to reconsider the whole analysis 
by first evaluating the force ( 40 ) in the zero-temperature limit at fixed 
d and then let d ^ co. In other words, the limits T — > and d ^ 00 
are not permutable, and the issue is about obtaining a simultaneous 
control of the force jointly for T near zero and d large. This will be 
the subject of forthcoming work. 



A Capacitor force 



Given that in the electrostatic part of the total force (36), only the 
monopolar part piP2dxiVAB of the loop Coulomb force dx^V^^ con- 
tributes (see Appendix |Bj), the electrostatic capacitor force fca.p{d) 
is 

/cap(d) = jdljd2jdye^,e^,pip2dx,VAB{l,2,y)pAil)pB{2). (86) 

The loop densities pa{^) and pb{2) are independent of y by space 
homogeneity in the plates' directions and Jdy dxj^VAB{^,2,y) = 27r 
(set q = in (41)). The remaining integrals factorize, and yield the 



particle densities pAixij'ji) and pb{x2,12) in plate A and B by means 
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of the identity 



p=i 



p{r,^) = yp D{X) p{C) (87) 



(see [37! Appendix D]). Introducing the charge density c(x) = Ylj e^/9(x, 7) 
in plate A and B, fca,p{d) is thus given by Formula ( pQ] ). Note that 
the charge densities ca{x) and cs(x) are still subject to the mutual 
interaction between the slabs, thus depend on the separation d. 
The magnetic part of the capacitor force is 

We show hereafter that 

jdy d,,WtB{l,2,y) = 1^ e*'=i(-i-2-rf) w-(xi, X2, ^i, fc = 0) 

(89) 

decays faster than any inverse power of (xi — X2 — d) as d ^ 00. This 
ensures that the decay of f^p{d) with the plates' separation has no 
power-law tail in view of 



The dipolar decay of >V™(1,2) (21) at large distance, responsi- 



ble for the power-law estimates (|81), is generated by the nonana- 
lyticity k^k^/K'^ due to 6^^^{k) in the transverse Coulomb potential 
4:TT5j];^{K) / k"^ {k = {ki,k)). However, setting k = eliminates the ki 
dependency in the transverse Kronecker symbol: 

.•-(fe.fc^Oj.^- iff-"^'- (90) 
10 if/i = i/ = l. 

Any nonanalyticity is thus removed in W^^{xi}X2, ki, k = 0) around 



A;i=0, ensuring the fast decay of (89). Indeed, in W'^{xi,X2,ki,k 







0), one is left with the stochastic integrals 

dri(si) • rdV2(s2)^^^^^Q(A;i,s~i-52)e*'=i^-i^i(^i)e-*'=i^-2^2(^2). 
Jo K 

Both g^{ki) and Q{ki 

, si— S2) are analytic functions of ki expandable 
as 1 -|- 0{kf). One sees by expanding the integrant around fci = 
that the only singular terms are functions of only si or S2- Their 
stochastic integration identically vanishes by Ito's lemma, stating that 
/(fdX(.)^O. 
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B Electrostatic force 



In Formula (32), the Casimir force has an electrostatic part /'^{d) due 
to OxVXb ^ magnetic part f^{d) resulting from differentiating the 
magnetic potential. One could write the average electrostatic force be- 
tween the two slabs directly by summing the Coulomb forces between 
the point charges: 



nd) 



lim 



lim [dri [dr2 e^^ey2dxVAB{ri,r2) 
.^oo L Jj^ ^ "—f 



71 72 



(91) 



where VABifa^fb) is the Coulomb potential (11), p\ is the particle 
density correlation function, and the same notation (31), (34) trans- 
lating positions from slab to slab B is used. Going to the phase 
space of loops by means of the identity 

oo oo « 

p!llL(^i>7i;r-2,72)= 5] Y^pm b{Xi) B{X2)pfs,Lih2) (92) 
(see [371 Appendix D] ) yields 

(93) 



In Formula (93), multipolar contributions of the Coulomb force are 
not present, in contrast to the electrostatic part of (33). The strict 
equivalence of these formulae relies on an invariance property regard- 
ing the choice of a reference point for a loop's position. Clearly, the 
loop 

£M = (rM,^,p,xM(-)) with xM(s) = X(s + u) -X(u) (94) 

describes the same path as the loop C = (r,7,p, X(-)) but has its 
origin r^""] = r + \^X(u) shifted by the vector A^X(n) (the time 
parameter is shifted by u). Such a shift does not affect the loop 
density: 

p(£M) = p(£) Vn, (95) 
whereas the two-loop correlation function satisfies 



(2),^ [mi] ^[U2h 

Pl y'-i ■''-2 I 



p^l\Cl,C2) if Ul - U2 G 



(96) 
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(see below). In the electrostatic part of (33 1, 



2) = [dsi ris25{si-S2) d^,VAB{r^t'\r^2'^) (97) 

JO JO 

(see ([20))) and one can replace l(^1) ^2) by j;^(£5*^^, £3 ^') (at 
fixed pi,P25 si, S2) because of the equal-time constraint in (97) that 
forces si — S2 to be integer. Performing first the changes of variable 



ri 



r2 



and then 



^i(-) 



(98) 



one obtains the electrostatic force 

lim -^/dl /d2e^,e^2 / dsi / ds25(si -s^) <9xi'yAB('^i, ■'^2)^^^^ 2)- 
L^ooL J A Jb Jo Jo 

(99) 



Indeed, the Jacobian of the transformations (98) is equal to 1: the 



random process X'-'^'{-) is still Gaussian with unit normalization, zero 



mean, and same covariance (17) as X(-), so that the Gaussian measure 
is unchanged: D(xM) = D(X) (see |371 Lemma 1] or [Ml Lemma 2]). 



Formula ( 99 ) eventually reduces to (91) since J^^dsi /(^Ms2'5(si — S2) 
Pi P2- 



The properties (95) and (96) both follow from the fact that the 



loop's self-energy in the activity (28) is invariant under a shift of origin 



and the loop pairwise interaction V{Ci,Cj) (19) is invariant when 



the loops Ci and Cj have their origin shifted to C^ ^ 



and Cf'^ with 

ui — U2 G The restriction ui — U2 € 7^ is the manifestation of the 
Feynman-Kac equal-time constraint in and of the quantum nature 
of the photon field in W™, occurring through the function Q{k, si — 
S2) (22): yV™ is unchanged because Q{k,{si + ui) — {s2 + U2)) = 
Q{k, Si — S2) when ui — U2 S Z by periodicity of the function s 1-^ 



C Screening of the resummed interac- 
tion ^ 

The classical Debye-Hiickel potential ^'^^^'^'^{xi,X2,y) for slab geome- 
try has been extensively studied in [26] . Because of the wall constraint 
on the screening clouds, this potential does not decay exponentially 
fast as would be the case in the bulk, but still has an integrable tail 
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~ along the wallj^ This imphes that its transverse Fourier trans- 
form <^'^'^^^^{xi,X2,k) is finite at fc = 0, see Formula (A. 11) of 
In the sequel, we infer that the same property remains true for the 



screened potential between loops $(l,2,fc) defined by (45). The only 
difference between the Coulomb potential v{ri — r2) for point charges 
and y^(l,2) is the extension of the loops that generates additional 
multipole interactions. To disentangle the monopole interaction from 
the multipole contributions, we proceed with the same method as in 
Sec. V.B.2 of [31] and only sketch the arguments. Introducing the 
multipole operator 

_^ _ / ^ Ijr^H^ .^^^2, (100) 

■^0 1=1 

the loop interaction is decomposed into its charge-charge (cc), charge- 
multipole (cm, mc) or multipole-multipole (mm) components: V^^ = 
+ + where 

Kc(l, 2) = pip2 v{ri - ra), Km(l, 2) = pi M2 v{ri - ra), 

¥^{-^,2) = M1P2 v{ri-r2), V^^{1,2) = Mi M2 v{ri-r2). 

(101) 

Summing the chains of V^^ to form <I> amounts to summing all pos- 
sible chains with bonds V^l, Vq-L^ ^mo ^-^d V^^. Summing first the 
pure V^^ chains builds the classical Debye-Hiickel potential $'=1'^**** with 
screening length k~^{x) = [47r/3 ^l-y I]p / D(X)p2 6^/3(2;, x)]"^''^j^ One 
is then left with the screened bonds 

i^cc = -/?e^ie^2 piP2 Fern = -/Je^.e^^ pi M2 

Frac = -Pe^,e^,MiP2 Fmm = -pe^,e^,MiM2 

(102) 

Finally, $ is built from convolution chains of these screened bonds 
subject to excluded convolution rules with respect to Fee. One sees 
that 

fP2 p 

Fe^(l,2) = -Pe^,e^,pi ds2 <I>=i^^^(ri, rs + A^,X2(s2)) - $=^'^^^(ri, rs) 

Jo L 

(103) 

is integrable in the y direction and so has a finite transverse Fourier 
transform at fc = 0. The same holds for the other screened bonds and 



^This was noticed long ago in Ref. see also [311 Sec. III.C.2] 
^This screening length reduces to the classical expression [47r/3^^ e^p(a;, 7)]^^/^ when 
exchange effects are disregarded [3TJ Sec. V.B.2]. 



30 



their chain convolutions, hence the result (47). These considerations 
apply to the various types of screened potentials considered in this 
work, e.g., the screened potential of the single plate systems 
and the potentials ^aa, ^bb occurring in Subsection 6.1. 



D Dipole potential Wab at order 0(d ^ 

From (Its]) and ([79]), the asymptotic form of Wa_b( 1,2, ^) is 



WAB(l,2,p 
1 ^^71^^72 



q . d~>oo 1 f dq 



n 

d y 27r * 

Pi 



P2 



dXf(.i)/ dX2^(s2) 



— ( 

27r 



+ q ■ Yiisi)] [qiX2{s2) + q ■ 1^2(^2) 



x=l 



-iqix 



ql + q^""" 



^%{qi.(i) 



x=l 



1 A^2 



Pl /•P2 
dXf(si)/ dX2^(s2 

Jo 



[iXiisi)i^ + q ■ Yiisi)] [iX2{s2)i: + q ■ Y2{s2)y{x, q) 



x=l 



(104) 



where 
v'"'{x,q) 



2-K qf + q 



TT 



-q\x\ 



X < 



^fiu + q\x\, 

iq^x, 



1^ = ^1=1, 

^i^l,u = l, (105) 



2S 



is the partial Fourier transform of the transverse Coulomb potential. 
The final result (Formula (5.88) in [40^) is obtained by working out 



the derivatives in (104). 
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